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A THREE-DIMENSIONAL ANALYSIS OF A TANGENTIAL YO-YO 


DESPIN DEVICE ON A ROTATING BODY 

By Robert L. Collins, Jr. 

Langley Research Center 

SUMMARY 

The problem of despinning a rotating rigid body by the release of yo-yo despin 
weights is considered. Although studies have been made previously for the two- 
dimensional case of a fixed spin axis, there has been little effort to incorporate an ini- 
tial coning angle which, in general, leads to a three-dimensional problem. This study 
considers the effect of an initial coning motion on the parameters which influence the 
design of a yo-yo despin system. The differential equations of motion of the system are 
derived for a symmetric rigid body and these are integrated numerically to provide solu- 
tions for several specific cases. The results show that a two-dimensional study is suf- 
ficient to predict such parameters as cable length, cable tension, and time required to 
despin, as long as the initial coning angle is less than 10° to 30°, and the masses of the 
despin weights are not too small. The parameters plotted in this study give an indication 
of the effects of initial coning angle on the design parameters of interest for what is 
called tangential motion. 


INTRODUCTION 

A typical problem area in the design of earth orbiting satellites is the reduction of 
the angular velocity which may occur after injection into orbit. One efficient technique of 
despinning a satellite is by the use of a so-called yo-yo despin device. Basically, this 
device is a simple arrangement of small masses attached to cables which are wound about 
the periphery of the spinning body in a plane perpendicular to the spin axis. When it is 
desired to reduce the rotational rate of the body, the despin weights are released and they 
move away from the body constrained only by the cables which unwrap from the body. 

The tension produced in the cables causes a moment which reduces the spin rate of the 
body. Essentially, the angular rotation of the body is reduced by transferring a desired 
quantity of angular momentum from the body to the despin weights. When sufficient 
angular momentum is transferred to obtain the desired reduced body spin rate, the cables 
are released from the body which then rotates in space independent of the yo-yo system. 



In this study it is assumed that the cables are initially wrapped in a circle about 
the body and therefore, during the unwinding phase of the motion, the cables will remain 
tangent to this ’’wrapping circle.” This phase of the motion has been referred to as 
’’Phase I” motion in reference 2; however, the terminology ’’tangential phase motion” is 
used in this paper. Many previous two-dimensional studies have also included a second 
phase to the motion in which the cables are assumed to have a fixed length and are 
attached to the body by a hinge at the end opposite the despin weights. Once the cables 
have unwrapped to their full length they swing away from the tangent until they are at 
some specified angle with the tangent, usually 90°, and then they are released. This 
phase of the motion has been referred to as ’’Phase II” motion. This paper does not con- 
sider despin cases which include this ’’nontangential” or Phase II motion. 

A particularly simple analysis of the motion of the yo-yo system may be made when 
there is no initial coning motion of the system. This motion is possible if the initial 
angular rotation vector of the system lies along the axis of symmetry. In this case the 
initial angular momentum also lies along the spin axis and if there are no disturbances 
on the despin weights or the body, the weights will move away from the body in a plane 
always perpendicular to the symmetry axis. The angular momentum vector, the angular 
rotation vector, and the axis of symmetry will then remain in a fixed direction in space 
and the resulting analysis can be treated as two dimensional. Thus, the two-dimensional 
motion requires that the movement of the despin weights and cables take place in a plane 
fixed in space and containing the initial circle of winding. 

A study of the tangential phase of the two-dimensional motion yields exact analytical 
expressions for the variables of interest. The equations of motion for the fixed or non- 
tangential phase in the two-dimensional case have not been solved as in the tangential 
phase but some results can be obtained through the conservation of energy and angular 
momentum principles. Studies of this type are presented in references 1 to 4. Refer- 
ence 2 gives a very concise treatment and also discusses the correction necessary to 
account, approximately, for the mass of the cables. Further studies of the two- 
dimensional yo-yo are found in reference 5 where it is noted that variations from the 
nominal design conditions may be partially nullified by using elastic cables. An analysis 
of this ’’stretch yo-yo” and some experimental data are found in references 6 and 7. 

An approximate solution for a three-dimensional problem is found in reference 8, 
but so many approximations are made on the configuration and kinematics of the system 
that it is not applicable for general studies of yo-yo despin systems. Reference 11, how- 
ever, gives a comprehensive study of the three-dimensional problem considering both 
phases of the motion but does not give results for complete despin by tangential phase 
motion only. 
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This paper presents a three-dimensional analysis of the parameters of motion 
which are of interest in the design and study of a yo-yo despin system. The assumptions 
made on the system are similar to those of the two-dimensional tangential phase studies 
with certain additional hypotheses which keep this analysis from becoming too complex. 
The equations of motion for the system are written by applying the Euler equations for 
rigid body motion to the rigid body and Newton's second law to a despin weight. This 
technique was used since it is directly applicable to this system while the Lagrangian or 
the Hamiltonian techniques require special attention and cannot be directly applied by the 
usual methods. The Lagrangian methods are invalid, in their usual formulation, for this 
problem because certain assumptions have been made about the configuration of the sys- 
tem which introduces constraining forces that do work during a virtual displacement of 
the generalized coordinates. This procedure is contrary to the basic postulates of ana- 
lytical mechanics and therefore the equations usually derived from this approach cannot 
be used without modification. In this vein it is also of interest to note that for the system 
studied herein, the total system kinetic energy decreases throughout the motion whenever 
the despin cables and weights are out of the initial plane of their winding. This decrease 
may be verified by making a detailed examination of the results presented and computing 
the system kinetic energy at several instants during the motion. Further discussion of 
this system may be found in reference 11. 

The equations of motion were integrated on digital computer and the solutions were 
checked by three separate means. First, the results for the two-dimensional theory 
previously mentioned agree with these solutions for the two-dimensional cases; second, 
the components of the system angular momentum vector were checked in fixed space and 
were found to remain constant, and third, an expression for the sum of the time rate of 
change of the system kinetic energy and the time rate of change of the work done by the 
system was found to be essentially zero as it should be. These checks are very con- 
clusive as to the validity of the results presented. 

SYMBOLS 


a radius of cable windings 

A acceleration vector of despin weight 

A nondimensional magnitude of acceleration of despin weight 

Aj,A 2 ,A 3 components of acceleration vector X along Xj, x 2 , and x 3 axes 
a^i elements of matrix from equations of motion (eq. (9)) 
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bi,b 2 , b 3 

bi,b 2 ,b 3 

Ci,c 2 ,c 3 

c 

E 1> E 2’ E 3 


i,j 




I 1’ I 2 ,I 3 


K 

l 

m 

M 


m b 


r 

t 

T 

T 


nondimensional form of elements a^j (eq. (12)) 
elements of inverse matrix of a^ 

auxiliary symbol for terms in equations of motion (eq. (9)) 
nondimensional form of bj,b 2 ,b 3 (eq. (12)) 

Coriolis and centrifugal acceleration terms (eq. (A13)) 
conversion factor (eq. (19)) 

Euler acceleration terms (eq. (A14)) 
indices 

unit vectors along xj, x 2 , and x 3 axes 
inertia ratio of body, I3/I1 

mass moments of inertia about xj, x 2 , and x 3 axes 

2m 

inertia factor for despin weights, 

h 

length of cable unwound 

mass of a single despin weight 

moment vector acting on body 

mass of body without weights 

distance from mass center to despin weight 

time 

tension in cable 

vector tension acting on body, directed toward weight 


4 





T 


nondimensional tension in cable 


V nondimensional velocity of despin weight 

components of velocity of despin weight along Xj, Xg, and X 3 axes 
(eq. (A7)) 

xj quantities representing /3, a, and T 

coordinate axes with origin fixed at mass center, directed along principal 
axes of body 

y v y 2 > • • *y 7 variables introduced to reduce second-order equations (eq. (14)) 
of cable out-of-plane angle (sketch ( 1 )) 

0 unwind angle (sketch ( 1 )) 

r coning angle (between k and angular momentum vector) 

6 spinning angle (between k and u>) 

nondimensional radius vector r/a from mass center to despin weight 

F 1 , r 2 , 5^3 components of r along x^, x 2 , and xg axes 

partial derivatives of r^r^fg with respect to 0 

r 1 ,r 0 ,r Q partial derivatives of r^,r 0 .v o with respect to a 
l a 6 ot ltd 

t nondimensional time 

Z) angular velocity vector 

a; nondimensional angular velocity vector 

o> nondimensional angular velocity vector (eq. (A 6 )) 

ojjjWg.Wg components of angular velocity along x^ Xg, and Xg axes 
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nondimensional form of 

w^2 nondimensional cross-axis spin rate, 

Dots over symbols indicate time derivatives; primes denote dimensionless time 
derivatives. A caret ~ over a symbol denotes a vector quantity and a tilde ~ over a 
symbol denotes nondimensionalization. Quantities subscripted with o indicate that the 
initial time (t = t 0 ) is being referred to. 

ANALYSIS 

Derivation of Equations of Motion 

The despin system which is analyzed in this report consists of a symmetric rigid 
body and two small masses or despin weights as they will be called. Each mass is 
coupled to the rigid body by a perfectly flexible, inextensible, and massless cable which 
is initially wound about the body. These cables are wrapped so that the despin weights 
are diametrically opposite each other through the mass center of the rigid body at the 
initial instant. The cables are assumed to be wrapped on the periphery of a circle which 
has its center at the mass center of the rigid body and is perpendicular to the symmetry 
axis of the rigid body. When the despin weights are released, they move outward and are 
constrained only by the cables which connect them to the body. The force of constraint 
acting on a weight is the tension in the cable and this tension manifests itself in a torque 
which tends to retard the spin of the rigid body. Therefore, there is a strong interaction 
or coupling between the motion of the rigid body and that of the weights. 

The body- axes coordinate system and some of the quantities used to describe the 
system configuration are found in sketch (1). A right-handed orthogonal coordinate sys- 
tem Xj,X 2 ,Xg is fixed in the rigid body with its origin at the mass center of the body, 
and the xg-axis lying along the symmetry axis of the body. Because of the symmetry of 
the body, any choice of orientation of the x^^ axes produces a set of principal axes. 
The choice for the x^-axis, however, will be made so that it is alined with a despin weight 
when that weight and cable arrangement is completely wound on the winding circle. The 
second despin weight then lies initially along the negative x^-axis. The xg-axis completes 
the right-handed triad. 

Since the despin weights are initially released with positions and velocities sym- 
metric through the origin, they will remain symmetric through the origin throughout the 
motion. Because of this symmetry it will be necessary to analyze only one despin weight 
and the weight initially along the positive xj-axis is used. 
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Sketch (1) 

The position of this weight is specified in terms of the body-axis coordinate system. 
As the cable unwinds, the point of contact where the cable leaves the body is defined by 
the vector a which is rotated an angle p from the xj-axis. The unwound cable is 
represented by the vector 1 in sketch (1) and must lie in a plane tangent to the cable 
winding circle and parallel to the x 3 -axis. The angle a is the angular measure of the 
vector 1 from the plane of the x 1( x 2 axes. It should be noted that, as defined here, a 
positive a corresponds to the despin weight being "below" the xj,x 2 plane in the sense 
of positive x 3 being "above" this plane. The angle p is referred to as the unwind 
angle and a as the out-of-plane angle. 

The length of cable unwound is given by the relation 

l = a/3 (1) 

where a is the radius of the cable winding circle and where l is the magnitude of the 
vector 1. With the constraint relation (1), the position vector r may be written with 
components along the body-axis coordinate system entirely in terms of the angles a 
and p. These components are seen from the geometry of sketch (1) to be 

rj = a(cos p + p sin p cos cf) 

r2 = a(sin p - p cos p cos a)\ (2) 

r 3 = -a/3 sin a 
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Because of the symmetry of the system, the mass center of the rigid body and the 
mass center of the system occupy the same position in inertial space. Therefore, the 
motion of the body is essentially that of a symmetric rigid body rotating about a fixed 
point at the mass center and being acted upon by a moment which is varying with time. 
This moment is due to the cable tension which acts along the vector i. The angular 
velocity vector u) of the body is denoted by the components u>j, u> 2 , and wg along 
the body axes. The force acting on the body due to the tension T in the cable is denoted 
by the vector T and the moment due to this tension is M. This moment is generated 
by the two cables and may be expressed as 

M = 2a x T 

The force vector due to the tension T and the vector a are seen from the geometry of 
sketch (1) to be 

T = T(sin /3 cos ai - cos /3 cos aj - sin ark) (3) 

a = a(cos /3i + sin /3j) (4) 


Where 1, j, and k are the usual unit vectors along the Xj, X2, and Xg body axes, 
respectively. The moment vector then becomes 

M = -2aT(sin or sin /3i - sin or cos /3j + cos ok) (5) 


The equations of motion of the body are obtained directly by use of the Euler equa- 
tions for rigid body motion about a fixed point. These equations may be found in refer- 
ence 10, and for this symmetric rigid body become 


I ld , l - (ij - l3)^2 a) 3 

*1^2 - ( ! 3 - I l) w l a, 3 
130)3 = -2aT cos a 


= -2aT sin a sin 
= 2aT sin a cos 



( 6 ) 


where I2 is replaced by its equivalent value Ij due to symmetry. 

Equations (6) are three differential equations for six unknown quantities a>j, 
cog, ol, /3, and T. Therefore, it is necessary to obtain three other expressions in order 
to specify the system completely. These expressions are obtained by applying Newton's 
second law to a single despin weight. 

The despin weight specified by the vector r, which has components along x^, X2, 
and X3 given by equations (2), has a force acting on it caused by the cable tension 
represented by -T. Therefore, 
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(7) 


-T = m 


d 2 r 

dt 2 


In this report the notation d/dt, d 2 /dt 2 refers to derivatives with respect to time in an 
inertial sense and (*), O') refers to derivatives taken in the Xj,X 2 ,Xg body-axis sys- 

tem. With this definition the absolute velocity and acceleration vector components along 
the body axes become 


dr - 

37 - = r + co x r 
dt 


, 2 ~ 
d r 

dt 2 




= ? + wXr + 2o>Xr + a)x(a)Xr) 


( 8 ) 


j 


Equations (2), (3), (7), and ( 8 ) may be manipulated to provide three second-order 
differential equations involving the quantities u>j, 0 ) 3 , Wg, a, p, and T and their 
derivatives. These equations will be independent of the three equations derived for the 
rigid body (eqs. ( 6 )) and therefore a complete set of six differential equations and six 
unknowns is established. 

The manipulations and differentiations required to obtain the desired relations are 
given in appendix A where it is shown that equation (7) may be written in the form 

3 

£ a ijXj = bj (i = 1, 2, and 3) (9) 

j=l 

The notation Xj which will be used again, signifies that the components of Xj are the 
quantities p, a, and T in that order for j = 1, 2, 3, and where the coefficients a^j 
are given functions of a and p only and the components bj are functions of wj, w 2 > 
o> 3 , a, p, &, and p. 

Equations ( 6 ) and (9) are the governing equations for this system, which can be inte- 
grated to yield o)j, cl> 2 , o> 3 , at, p, and T. These equations are nonlinear and it is 
very doubtful that an exact closed form solution can be obtained; therefore, integration of 
the equations by some numerical procedure is necessary. 


Nondimensionalization 

In this section the parameters which appear in equations ( 6 ) and (9) are put into a 
convenient nondimensional form for computation and for presentation of the results. It 
may be noted that equations ( 6 ) and (9) as they stand contain the variables coj, ajg, ^ 3 , 
a, p, T, and t and constants Ij, I 3 , m, a and that after nondimensionalizing, the 
variables are Sj, £ 3 , ^ 3 > T, and t whereas the constants are I and K. 

Note that the number of variables is the same and the constants are reduced to two. 
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When nondimensionalized in the manner presented herein, the variables become 
independent of the constants I and K for the tangential two-dimensional case. This 
form then has the advantage that one curve, in the family of curves to be presented, cor- 
responding to the two-dimensional or zero coning case will be the same for any of the 
values of the I and K which are to be chosen for study. This choice of nondimension- 
alization therefore results in more convenient and uniform plots. For these reasons, the 
following nondimensional quantities are defined. 


I 3 

'-il 




K = 


2 ma 2 


w, 


wi = 


co 


3o 


CO, 


w 2 - W 


3o 


CO, 


co Q = — 
3 co 


/3 = 


3o 


jm 

V K 


a = 


a 


1 + I 


K 


■Rt 


2aT 


3 W 3o 


co 


r = 


3o 


1 + i 


K 


(10) 


Note that 


K " 2 ma 2 


The subscript o in implies that the angular velocity at the initial time is being 

considered. 

The equations of motion can then be rewritten in terms of these new parameters. 
If the derivatives with respect to the nondimensional time r are denoted by primes, 
equations ( 6 ) become 
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( 11 ) 


s; = (i - 1 )\j- 


i ~ 


1 + ^ a) 2° L) 3 “ IT sin a sin Z 3 


= (I - D|/l +^2 1 5 3 + IT sin of cos /3 


Wg = -T cos a? 


where a and /3 are 


a = a)jl + ^ 

e = ?f71 

equations (9) rewritten in terms of these nondimensional parameters become 

3 

V a^.X. =b. (i = 1, 2, and 3) (12) 

j=l 

where again Xj corresponds to quantities /T , 5 ", and T and where a^ and bj 
are functionally 

bi = ^ t K;Z v Z 2 ,Z z ,a',p,a,p) 

The complete expressions for the a^ and b^ are found in appendix A along with a 
derivation of the equations and their nondimensionalization. Equations (11) and (12) are 
the equations of motion which will provide the information to be presented in this report. 

Integration of Equations of Motion 

Equations (12) are not directly amenable to the usual techniques of integration by 
digital computer since the left-hand side is a linear combination of the variable T and 
of the second derivatives of 5 and j3. Therefore, in order to put these equations into 
a more fundamental form, the components Xj are solved by use of Cramer's rule so 
that 

3 

Xi = V iijbj (i = 1, 2, and 3) (13) 

j=l 
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where the ajj are the elements of a matrix which is the inverse of the matrix with 
elements a^j. 

Equations (13) provide explicit relationships for the derivatives a" and jgf' and 
also for the variable T. In order to obtain a complete set of first-order differential 
equations in place of equations (11) and (13), a new set of variables are defined as follows: 


yi = s i 

y 2 = “2 

ys = s s 


y 4 = P > 

y s = 5 

y 6 = P 


(14) 


y 7 = aj 

With these new variables, equations (11) and (13) may be written as a set of algebraic and 
first-order differential equations: 


3 

9 * I r 3i E i 

i=l 

y i = Jl + - I)y 2 y 3 “ sin^l + i a sin^|l + i /3 

y 2 = l/ 1 + " l) y 1 y 3 + IT sin ]jl + i a cos^l + i f} 

y 3 = - fcos f+];* i 

y 4 = y 6 

y 5 = y 7 
3 

yg = J a^b. 

i=l 

3 

y? = I hA 

i=l J 


(15) 


12 



Equations (15) are now in a convenient form for numerical integration. However, 
at the initial instant the coefficient matrix of equations (12) is singular. Hence, special 
methods must be used to begin the integration and a simple technique is shown in appen- 
dix B which overcomes this difficulty. The only initial quantities needed to begin the 
integration are the constants I and K and the variables and ti>2 0 * other 

variables are obtained by setting the unwind angle /3 and the relative velocity r equal 
to zero at the initial instant. This result is also shown in appendix B. 


Descriptive Parameters of the Motion 


Thus far, it has been shown that the motion of the tangential yo-yo despin system 
analyzed in this report can be described in terms of the quantities cDj, cDg, 5, /3, 

the nondimensional tension T, and the nondimen sional time r. A particular motion 
results from a choice of the quantities I and K and the initial conditions c5j o and 
u>2o’ Two other quantities which are of interest are the initial coning angle r 0 and the 
cross-axis spin rate The initial coning angle is the angle between the angular 

momentum vector and the axis of symmetry of the system before the despin process 
takes place. The motion before deployment of the despin weights is characterized by 
this coning angle and is described in references 8 and 9. 


Since both the axis of symmetry and the angular momentum vector of the rigid body 
move in inertial space during the despin process, the instantaneous value of the coning 
angle is somewhat academic, but the initial value r o is important in the physical inter- 
pretation it provides. The coning angle is computed from 


r = 



(16) 


where ^12 * s cross-ax i s spin ra te an ^ I s the magnitude of the nondimensional 
angular rotation vector in the x^Xg plane or 


S 12 = (17) 

Another angle which has some academic importance and which is presented in the data is 
the spinning angle 6. This angle is the angle between the angular rotation vector and 
the body symmetry axis or 


6 = 


tan 



(18) 


From the relations (16) and (17) it may be seen that the initial values and u>2 0 

may be determined from r o and ^o/^lo' * n analysis these latter 
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quantities are used as initial conditions rather than 
the relations necessary to obtain and u> 2 0 from 


and £> 2 q - Appendix B gives 
r o and G lo/ G 2o- 


RESULTS AND DISCUSSION 


The equations of motion can be integrated once the constants I and K and the 
initial conditions ^o/^lo anc * ^o are specified. The integrations are ended when 
the symmetry axis spin rate is reduced to zero or if some other condition which is 

physically unrealistic occurs. For instance, the integrations end whenever a becomes 
90°, since this condition would most likely correspond to physical interference between 
the cable or despin weight and the rigid body. Another physical restriction on the vari- 
ables occurs when the unwind angle /3 reaches a maximum value and begins to decrease. 
This "rewinding” of the cable (that is, /3 decreasing) is not physically possible since it 
requires a constraint at the point of contact of the cable and winding circle (defined by 
a) to "tuck" the cable back into or along the winding circle. Since a mechanism to pro- 
vide such a constraint is not considered in this study, the event of /3 decreasing is not 
allowed. 


Three values of I were chosen for this investigation which would represent dif- 
ferent geometrical configurations for a rigid body of constant density. Sketch (2) shows 
the relative geometry of right circular cylinders and cones which have the inertia ratios: 
I = 0.1, 0.9, and 1.3 used in this study. 


Xr 



I = 0.1 



x 3 





Sketch (2) 

The inertia ratio K is an important design parameter since it contains both the 
despin mass m and the initial cable winding radius a. In order to obtain realistic 
values for this parameter, the right circular cone and right circular cylinder were again 
considered. If a constant mass density is assumed, it is possible from the equations of 
mechanics to derive the following relationships: 
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(Cylinder) 


K = 



K 


- 15 /m \ T 

= 4 Kj 


(Cone) 


where m^ is the mass of the rigid body and m is the mass of a single despin weight. 
Therefore, if I is known, the inertia ratio K is determined by the ratio m In 
this paper the value of the mass ratio was chosen to lie between the limits of 0.0005 to 
0.05 which correspond to a total mass for both despin weights of 0.1 percent to 10 percent 
of the body mass. The value of K from these considerations must be approximately 


0.0002 JKi 0.02 
0.0018 S K £ 0.18 
0.0026 S K £ 0.26 


(I = 0.1) 
(I = 0.9) 
(I = 1-3) 


The parameter ^o/^lo was l a ^ en as unity in every case of this study. This 
value was assumed to be a representative one and a change in this value should not affect 
the general trend of the results obtained here, especially for the smaller coning angles. 
Changing this ratio essentially varies the initial conditions a°, ft' Q , and T 0 . 

The initial coning angles studied are r Q = 0°, 10°, 20°, 30°, 40°, and 50°. The 
angle of 50° seems to be a reasonable upper limit for practical considerations since at 
this value the body has a larger component of angular momentum in the x^Xg plane 
than along the symmetry axis. It should be recalled that the initial condition r 0 = 0° 
implies that the ensuing motion is two dimensional. 

Figures 1, 2, and 3 are time histories of the quantities Sjg? £>3, T, T, 6, 5, 

and /3 for various values of r Q and for three specific sets of values for I and K. 
Figures 4, 5, and 6 are presented to give some indication of the effects of changes in the 
design quantity K on several important parameters of the motion for the three geome- 
tries specified by I = 0.1, 0.9, and 1.3. 


Time histories for a particular body despin weight system I = 0.1, K = 0.001 
are presented in figure 1. The most significant features are the nearly constant values 
of the cross-axis spin and small effect of coning angle on the spin rate, cable tension, and 
cable length. The design quantities ct>g, ¥, and /3 could have been predicted very 
accurately by the two-dimensional techniques given in reference 2. In figure 1(d) the 
coning angle and spinning angles are presented. At despin = 0), these angles are both 
90° which will always be the case when some residual tumbling rate or residual cross 
axis spin rate exists. Figure 1(e) reveals that the maximum value of the out-of-plane 
angle is approximately two- thirds of the initial coning angle. 
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In figure 2, time histories are presented for the same parameters as in figure 1 
but for the different geometrical case I = 0.9, K = 0.01. It is significant to note that 
there is considerable deviation from the two-dimensional case of all the variables plotted 
for initial coning angles of about 30° or larger. Also note that there is a decrease in the 
cross-axis spin rate during the despin which is usually desirable since the final 

value is the residual tumbling rate after despin. The out-of-plane angle (fig. 2(e)) has a 
maximum value which is approximately twice the initial coning angle. 

It was mentioned previously that whenever the angle (3 reached a (local) maximum 
during the despin process and then began to decrease, the integration was stopped. This 
event occurs on the T 0 = 50° case of figure 2. Figure 2(f) shows j3 reaching its 
maximum where its slope is zero. At this point the integration is stopped and this occur- 
rence is noted in figure 2(b) which shows that the symmetry axis spin rate has a value of 
about 0.075 at the end of the integration. 

Figure 3 is a set of time histories of the significant parameters of the motion for 
I = 1.3 and K = 0.015. The trends of the variations are similar to those of figure 2; 
however, the magnitudes of the variations are considerably larger. The maximum values 
of the angle a are, as in the previous case, about twice the initial coning angle. The 
case of (3 reaching a maximum value does not occur in these cases, although, as indi- 
cated in figure 3(f), the slope of the curve of 0 plotted against r for r Q = 50° and 
t = 0.5 is very nearly zero. When compared with the curves of figure 1, the variations 
of these curves from the two-dimensional case (r Q = 0°) are very large. 

The figures showing the nondimensional cable tension (figs. 1(c), 2(c), and 3(c)) 
have been nondimensionalized by the quantity Wg o and not by the total rotational velocity 
magnitude u> 0 . A false impression of the effect of coning angle on the cable tension may 
be given if these figures are not interpreted properly. If it is desired to obtain the change 
in cable tension because of a change in the coning angle for a fixed value of the magnitude 
of the angular rotation vector, a nondimensionalization based on the quantity o> 0 should 
be used. A simple conversion for altering the nondimensional form is 


where 


2a\ 



K 


l3^o 2 


CT 


C = 


1 

1 + (I tan r Q ) 2 


(19) 


If this conversion factor is used for the case presented in figure 3(c) and for t = 0.6, the 
value of nondimensional tension based on co 2 rather than oj 3o for r 0 = 50° is 0.7 as 
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compared with 2.35 for T. Therefore, when compared with the r 0 = 0 case, it is 
observed that for a given value of the total rotational magnitude u> 0 , the cable tension 
decreases with increasing initial coning angle. 

In figures 4, 5, and 6, four parameters of basic interest for the design of a despin 
system are presented as a function of the ratio K/l for various values of the initial 
coning angle. These curves represent the three basic geometries I = 0.1, 0.9 and 1.3. 
Values of T > /3 at the instant when is zero are presented along with the 

maximum value of the nondimensional cable tension. 

In figures 5 and 6 an absence of data for T 0 = 50° may be noticed. This lack is 
due to the occurrence of the rewinding phenomena (/3 = /3max) mentioned previously and 
shown in figure 2. Also, in figures 5 and 6, some of the curves are dotted and some 
dashed for reading clarity in the rather congested area at low values of K/l. 

The quantity a t Wg = 0 represents the residual tumbling rate of the system 

and its variation with T 0 and K/l are given in figures 4(a), 5(a), and 6(a) for the three 
values of I considered. For I = 0.1, this final cross-axis spin rate is independent of 
the inertia ratio parameter for all values of the initial coning angle up to 50°; however, 
for I = 0.9 and 1.3, the design problem is a little more sensitive, especially at the lower 
values of K. 

The nondimensional despin time is shown in figures 4(b), 5(b), and 6(b). This 
parameter can be important, for instance, when the despin of spinning atmospheric entry 
probes is being considered. It is seen again that for 1 = 0.1 the two-dimensional despin 
time (r Q = 0°) is a very good estimation of the three-dimensional despin time except 
when K/l is chosen very small and r 0 is large. For the larger values of I, the 
despin time would be grossly in error if the two-dimensional results were used where 
there was a significant amount of coning. 

The maximum nondimensional cable tension shown in figures 4(c), 5(c), and 6(c) is 
also very near the two-dimensional value for I = 0.1 except when K/l is very small, 
whereas at the larger values of I both initial coning angle and changes in K produce 
large variations in this parameter. 

The parameter /? at ci>g = 0 represents the nondimensional cable length required 
for despin and is plotted in figures 4(d), 5(d), and 6(d). The variations of this parameter 
for changes in T 0 and K are again small for I = 0.1 but grow larger for the larger 
values of I. 

In general, figures 4, 5, and 6 show that for a slender body (I = 0.1), the two- 
dimensional theory is probably sufficient for most engineering purposes unless very small 
despin weights or very small values of a are desired, whereas for more disk-like bodies 
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(I = 0.9 and 1.3), a more thorough analysis should be made, especially when large coning 
angles or small values of K are to be encountered. 


CONCLUDING REMARKS 

Integration of the equations which describe the three-dimensional motion of a 
spinning rigid body during the yo-yo despin process must be done numerically because no 
general analytical solutions to these equations have been derived. This paper presents 
the equations and their derivation and points out that special precaution must be taken at 
the beginning of the numerical integration procedure because of a singularity in the coef- 
ficient matrix of the differential equations at the initial instant of motion. A convenient 
method of nondimensionalization is also presented which allows the reduction of the 
parameters necessary for integration to be reduced to four: inertia ratio of body I, 

inertia factor for despin weights K, initial angular velocity ratio ( ^2o/ co lo’ and 
coning angle r 0 . An analysis of some particular cases shows that despin time, cable 
tension, and cable length may be computed with good accuracy by the two-dimensional 
approximations for some range of F 0 and K/l for all the values of the inertia ratio I 
studied. For I - 0.1, it is found that the two-dimensional approximations hold very well 
for all T 0 if the ratio K/l is above 0.02. For I = 0.9 and 1.3, it is found that the 
range of applicability of the two-dimensional results should be restricted to values of 
r 0 of about 10° to 20° depending upon the accuracy which is required. 

In every case studied herein, the cross- axis spin rate after despin is found to be 
less than or equal to its initial value, and therefore, it appears that the yo-yo device also 
has a damping effect on this component of the angular velocity vector. 

For this study of the tangential despin, the value of the angle a did not increase 
beyond the value of the initial coning angle for the three time histories presented and did 
not reach a value of 90° for any of the cases integrated. 

An important point should also be made here concerning the length of cable required 
to despin the body. It can be observed that the nondimensional cable length ft required 
to despin is very sensitive to both the coning angle and the ratio K/l. A slight misjudg- 
ment of this value could lead to serious difficulties if some initial coning exists. For 
instance, if the cable length is designed from two-dimensional considerations, it may be 
too long and may not unwind and release properly. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., September 7, 1966, 

124-07-02-35-23. 
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APPENDIX A 


DERIVATION OF THE EQUATIONS OF MOTION OF A DESPIN WEIGHT 

In the paper, it was pointed out that the equations of motion of a despin weight can 
be written as 


-T = m 


where 


r = r^i + rgj + r^k 


and Tj, r^, and r^ are given in equations (2) and T is found in equation (3). In 
this appendix the nondimensional forms of these parameters are used and, therefore, 
equation (12) of the main text is derived here rather than equation (9). 

The nondimensional radius vector r is formed as 

When derivatives are taken in terms of the nondimensional time r in place of t, the 
following substitutions may be made 


O-O’Sf 

dt 


or using the definitions of equations (10) 


(•) = ()' 


i + i 


(••) = ()" 


l + i 


Equation (7) may be rewritten as 


JdLt. =raa ‘>4— 

2af7I 1*1 
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~T - K d 2 ? , A5 . 

It is found convenient for computations to introduce another nondimensional form for the 
angular rotation vector o> 

c* = ^ s (A6: 

The derivatives of the vector r can be formed in the following manner: 

V = ^ = r + u)* x r 


V 1 " ? i + *2 r ~3 - S 3 ? 2 
V 2 = ? 2 + S 3 f l - S l f 3 
V 3 = r 3 + *1*2 - w 2 r l 


The relative derivatives of the components of r are found from 


r\ = ft S’ + F; a' 

1 X /3 X Q! 


(i = 1, 2, and 3) (A8) 


and from equations (2) of the main text 


*lp s Tf =sin P cosa -*t 


'2e = dp~ = *i- cos 0 cos « 


ro = -rx- = -sin a 

J /3 9/3 


1 ~ 


'l a *Tta =r 3 Sln < ! 


? 2„ * U = ' ? 3 cos f> 


To = -j3 cos a? 
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The "absolute" acceleration is then 


A l2^ A A ys 

t dr ~ ^ '‘ w 

A = — x = r + co x r + tn x r + u> x V 

dT 2 


(A10) 


The "relative” acceleration of r is needed and is 


r" = r.” = r. p ' + r. a" + r. ' p' + r. V (i = 1, 2, and 3) (All) 
1 x p 1 a 'p 'a 


where 


r = |cos p cos a - *2^' ~ ^ s * n ^ s ^ n 01 + ^2^' 

? 2 ^' = ^r + sin p cos otjP' + ^r ^ + cos p sin aja' 

rg^' = -(cos a)a' 

? 1 a = (% sin a + f 3 cos py + (% sin p) a ' 

% = ( F 3 sin 0 - % cos - (% cos 


r, ' = -(cos o?) (3' + ( p sin a)a' 
°a 


(A12) 


In equation (A10) the Coriolis and centripetal accelerations may be written together as 


Cj where 


and the Euler terms 


C 1 ~ _ w 3 r 2 + ^2^3 “ w 3^2 

C 2 = w 3 r l " w l r 3 + W 3 V 1 - W 1 V 3) 
C 3 = u) l r 2 _ Cl, 2 r l + S 1^2 " ^2^1 


(A13) 


Ei = 5|’r 3 - S>3'r£| 


E 2 = r ^ - 5 Jr 
E 3 = W 1 r 2 " w 2 r l 




(A14) 


If the Euler equations are used as given in equations (11) along with the definition in equa- 
tion (A6), the derivatives £>*' may be replaced in equation (A14) by the variables 53;, 

P, a, and T. Also if the main text of this paper is referred to for the components of T, 
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the derivatives and other expressions derived in this appendix may be manipulated to 
obtain an expression in the form of equation (12) or 

3 

Y ajjXj = q 
i=i 


(i - 1, 2, and 3) (A15) 




where corresponds to ft , a , and T and 


a ll = r D 


a 12 “ r l„ 


~ sin j3 cos a + Ir^ sin a cos ft + rg cos a 


a 21 = v 2r. 


a 22 " r 2 a 

a 23 = ~ cos ^ cos 01 ~ ^1 COS 01 + ^3 s * n a s * n ^ 


a 31 = r 3, 


a 32 " r 3 


a 


a 33 ~ ~ K S * n 01 ~ ^2 s * n 01 sin ^ ” ^1 s * n a cos 0 


b l = 


- C 1 - ? l/ " ? 1 + ? 3 (1 - I)53 1 S 3 


L Q! 


R 


bo - 


- * 


C 2 - r 2ff' ~ T 2f + r 3 (1 - I)w 2 w 3 


R 


-C 3 - ?g j3* - rg a ' - r^l - IJcDgSg - r^(l - I)<ij£g 


1 + 


K 


where a and j3 may be written in terms of « and /3 by use of equations (10) and go* 
may be rewritten by equation (A6). These expressions give the information required to 
integrate the governing equations (12). 
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INITIAL CONDITIONS 


The conditions which must be satisfied at the initial instant are: (1) the velocity of 

the despin weights relative to the body is initially zero and (2) the unwind angle j3 is 
also zero. The angle a at the initial instant is not prespecified by the choice of coor- 
dinates. However, when the conditions r Q = 0 and /3 0 = 0 are applied, it is found the 
a Q = 0 is a necessary initial condition. Furthermore, when these conditions are applied 
to the coefficient matrix in equations (12), all the elements a^j are found to be zero 
except agg which is the coefficient of the initial nondimensional cable tension T in 
the second equation. Equations (12) then provide the following expressions which must 
be valid at the initial instant: 


K - p 


2o) 


+ 1 


T = 
1 o 


"lo w 2o 


1 + 1 


K 


“0 = 


2 - I 

2£n 


CO 


lo 


(Bl) 


If it is noted that 

S 3o = 1 

it can be observed that if I, K, u5j 0 , and cDgo are known, the initial values of all the 
other quantities Sq, |3 0 , a Q , fi' Q , and T 0 are specified. The fact that only four quan- 
tities are necessary to specify the motion of this system leads to reductions in the amount 
of effort required to present the data. 


Equations (15) could now be integrated if it were not for the fact that the matrix of 
elements a^j is singular at the initial instant and therefore T, y^, and y^ cannot be 
computed initially. It is, therefore, necessary to find other expressions which will be 
valid at the initial instant. This procedure has already been done for T 0 as was shown 
in equations (Bl). The initial values of yg and y^, which correspond to /3 Q and « 0 , 
are found by differentiating equations (12) with respect to the nondimensional time r and 
applying the initial conditions which are given. When this is done, two relations valid at 
the initial instant are found 


— * T ~ 

a ll^o + a 12°'o = b i " a i* T r 


13 o 


I /s/1 . 1 __ 

a 31^o + a 32°o = b 3 ■ a 33 T o 


(B2) 
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After taking the appropriate derivatives and manipulating, the desired initial conditions 
are found to be 


(B3) 


which are the initial values of the quantities y g and y ^ in equations (15), and the 
values of and a5g Q may be found from equations (11) at the initial instant. 

These sets of initial conditions are used to start the integration and once started, equa- 
tions (15) are integrated by a standard fourth-order Runge-Kutta numerical integration 
technique. 

In place of presenting initial conditions in terms of u>j o and a>2o’ ^ was decided 
that the initial coning angle r Q and the ratio &2o/^lo woul<1 be more significant and, 
therefore, these quantities are the parameters which represent the initial conditions in 
the data presented. If r 0 and ^o/^lo are £* ven initial conditions, the quantities 
o>l 0 and uigo must be determined from them so that equations (15) may be integrated. 

By noting that 

^2o = ^2o 
5 lo _a; lo 

and by equations (16) and (17), the component of the nondimensional angular velocity vector 
along the Xj body axis at the initial instant is 


and the Xg component 






(B5) 


The e in equation (B4) is chosen as positive or negative unity so that the com_ 

ponent of the angular rotation vector is put into the appropriate quadrant in the x^ Xg 
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plane. For this study it is assumed without loss of generality that the Xj, Xg, and x^ 
coordinate system is chosen so that at the initial instant, the parameter C 2 Q is posi- 
tive. In this case e is determined by 

6 = Sgn fe) (B6) 
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